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Formalism for radiative correction (RC) calculation in exclusive pion electroproduction on the 
proton is presented. A FORTRAN code EXCLURAD is developed for the RC procedure. The 
numerical analysis is done in the kinematics of current Jefferson Lab experiments. 
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Understanding the electromagnetic transition ampli- 
tudes from the nucleon ground state to excited states 
provides valuable insights into the electromagnetic struc- 
ture of the nucleon. Exclusive pion electroproduction is 
one of the major sources to provide the most direct infor- 
mation about the spatial and spin structure of the excited 
states. With the development of high intensity and high 
duty-factor electron beam with high degree of polariza- 
tion, this field reaches a new level of quality. For the past 
several years, exclusive pion electroproduction has been 
the main subject of extensive studies at various acceler- 
ator laboratories such as MIT-Bates, ELSA, MAMI and 
Jefferson Lab. 

New measurements with the CLAS detector at Jeffer- 
son Lab/Hall B (lj are expected to greatly improve the 
systematic and statistical precision and cover a wide kine- 
matic range in four-momentum transfer Q 2 and invariant 
mass W, as well as the full angular range of the resonance 
decay into the nucleon-pion final state: 

e(ki)+p(p) — > e'(fc 2 ) +TT + (ph) +n(p u ), 
e{ki)+p{jp) — ► e'(k 2 )+p(p h )+Tr"( Pu ), (1) 

where Ph(Pu) denotes the momentum of the detected 
(undetected) hadron. Adequate calculation of radiative 
corrections (RC) becomes important in interpreting the 
measured observables such as unpolarized coincidence 
cross sections and polarization asymmetries. 

While solving the RC problem, people are commonly 
referred to the classical approach developed by Mo and 
Tsai g| and used for inclusive and elastic electron scat- 
tering for decades. However, this approach cannot be 
directly applied for exclusive pion electroproduction due 
to the following reasons. 

First, we now deal with exclusive electroproduction, 
where the hadron is detected in addition to the final elec- 
tron. It reduces the room for phase space allowed for the 
final radiated photon. The formulas of Mo and Tsai as 
well as any other inclusive formulas cannot be applied for 
this case without additional strong assumptions. 

Second, there are only contributions of two structure 
functions in the inclusive case. The exact formalism for 
the exclusive process requires consideration of four struc- 
ture functions for the unpolarized case with additional 
angular dependence associated with them. Note that 
transition to the case of two structure functions would 
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FIG. 1: Feynman diagrams contributing to the Born and the 
next-order electroproduction cross sections, a) Born process, 
b) and c) Bremsstrahlung, d) Vertex correction, and e) Vac- 
uum polarization. The momentum ph(pu) is assigned to the 
detected (undetected) hadron. 



not be possible even within realistic approximations. Mo 
and Tsai's approach predicts neither RC to polarization 
asymmetries, nor dependencies on the outgoing hadron 
angles. 

The third reason is a known shortcoming of the Mo 
and Tsai's approach, namely, the dependence on an un- 
physical parameter splitting soft and hard regions of the 
phase space of radiated photon in order to cancel the 
infrared divergence. 

In our approach, which is based on a covariant pro- 
cedure of infrared divergence cancellation proposed by 
Bardin and Shumeiko in Ref.||, such an unphysical pa- 
rameter is not required. Previously, this approach was 
applied for the calculation of RC for inclusive [ I , [| [| , 
semi-inclusive J?], |j and exclusive diffractive [£] reac- 
tions. Recent reviews of the approach, higher order 
effects and calculation for specific experiments can be 
found in papers 0, |ll|, |r3 . Based on these results, a 
Fortran code POLRAD [Of for RC calculation in po- 
larized inclusive and semi-inclusive processes was devel- 
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oped. Besides, some specific tasks such as Monte Carlo 
generator RADGEN [lj], MC approach to diffractive 
vector meson electroproduction [p~5j , RC to spin-density 
matrix elements in exclusive vector meson production 
[M, quasielastic tail for polarized He-3 target were 
solved. Recently Bardin-Shumeiko approach was applied 
to the measurements of elasticpolarized electron-proton 
scattering at Jefferson Lab Jl8|, |l9). A comprehensive 
analysis of results obtained in Ref. Q and [|| was made 
inRefs. @,|2^]. 

Note that there are also other approaches for calcula- 
tion of RC in electroproduction processes. For example, 
the results of Ref. are actively used for HERA ex- 
periments, while the approach developed in papers J22] 
is applied to specific measurements at Jefferson Lab [E3f. 

The Feynman diagrams needed to calculate RC are 
presented in Fig.|l|. They include QED processes of radi- 
ation of an unobserved real photon, vacuum polarization 
and lepton-photon vertex corrections. These processes 
give the largest contribution due to a large logarithmic 
term m(Q 2 /m 2 ). They can be calculated exactly from 
QED, and uncertainties of such a calculation are only 
due to the fits and data used for the hadronic structure 
functions. These uncertainties are demonstrated in the 
present article. Additional mechanisms (box-type dia- 
grams, emission by hadrons) are smaller by about an 
order of magnitude and they contain considerable theo- 
retical uncertainties. Most recent studies of two-photon 
exchange effects for elastic ep-scattering were reported 
in Refs.j24|, [25). Previous RC calculations |^(| for inclu- 
sive polarized deep-inelastic scattering included above- 
mentioned additional mechanisms. Generalization to the 
case of meson electroproduction will be subject to a sep- 
arate study. 

The paper is organized as follows. We introduce kine- 
matics and definitions (Section II), derive a cross section 
of the radiative process (Section III), solve the problem of 
infrared divergence (Section IV), obtain RC in the lead- 
ing log approximation (Section V), verify the relation be- 
tween exclusive and inclusive RC (Section VI), perform 
numerical analysis (Section VII), and summarize the re- 
sults in Section VIII. 




FIG. 2: Definition of momenta and angles in center-of-mass 
frame. 



Ui = S - Q , u 2 = X + Q , Ai,2 = u 12 
S P:X — S ± A, As = S 



4m 2 W 2 , 



2 4m 2 M 2 



w z = s x 



M 2 ,X q = Si + 4M 2 Q 2 , 



X = Q 2 u lU2 - Q 4 W - m 2 A„ 



(2) 



where m(M) is the lepton (proton) mass. 

We use the cm. system of virtual photon (q — k\ — 
k 2 ) and initial nucleoli. The axis OZ is chosen along 
q (see Fig. 2). The energies and angles (Fig. 2) in the 
selected frame can be expressed in terms of invariants 
(Q) as follows: 

p _ m,a p _ S X ~ 2Q 2 _ S X + 2M 2 
*i,a 2W q ~ 2W ' p ~ 2W ' 



P1.2 = 



-, Pp=Pq 



COS Pi 2 = 



2W ' 1H 2W 1 

ui, 2 (S x - 2Q 2 ) ± 2Q 2 W 2 



II. KINEMATICS AND BORN PROCESS 

At the Born level (Fig. pi) , the cross section of the pro- 
cesses (Q) is described by four kinematic variables. Fol- 
lowing tradition, we choose them as squared virtual pho- 
ton momentum Q 2 , invariant mass of initial proton and 
the virtual photon W, and detected pion (or proton) an- 
gles Oh and 4>h in the center-of-mass of the final hadrons. 
The obtained formulas are equally valid for both electron 
and muon scattering. 

We use the following Lorentz invariants defined from 
leptonic 4-momenta: 



S = 2k lP , X = 2k 2 p, Q 2 



fe) 2 , 



sin 0i,2 



2WV\ 



(3) 



Here subscripts 1,2, q,p denote the initial lepton, final 
lepton, virtual photon and initial proton, respectively. 

All the introduced kinematic variables have the same 
definition for both Born and radiative kinematics. But 
the situation is different for the final hadrons. The rea- 
son is that the energy of the observed hadron is not fixed 
by measurements of the chosen kinematic variables: Q 2 , 
W 2 , Oh <j>h. As a result, it is different for these cases. In 
the Born case the hadron energy can be defined via con- 
servation laws, but for the radiative process it depends on 
the unknown energy of the radiated photon. In order to 
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distinguish between the variables, we use the superscript 
'0' for the quantities calculated for Born kinematics: 



W 2 



Ph 



2W ' rn 2W ' 
\° w = (W 2 -m 2 ~m 2 u ) 2 ~4m 2 W 2 . 



(4) 



All invariants defined via the measured 4-momentum of 
hadron ph are also different for Born and radiative cases: 



V£ 2 = 2k h2 p h = 2{E ia E° h - pi, 2 p°(cos h cos h2 
sin Oh sin 0i )2 cos <j> h )), 



2 PPh = W 2 + m 2 h -m 2 u + V 2 ° - V, . 



(5) 



The last expression follows from 4-momentum conserva- 
tion. 

The Born cross section of exclusive pion electroproduc- 
tion in terms of the introduced variables reads: 



d(J 



M 2 



dk 2 dp u dp h 



2S(2tt) 5 2E 2 2E u 2E h 
(Ana) 2 L%W^ 



2(4tt) 4 5 2 Q 4 



S(A-p u ) 



Xw dQ 2 dW 2 dQ h ,(6) 



W 2 



where A is the total 4-momentum of the undetected par- 
ticles, A = p + k\ — k 2 — ph. For the Born process 
(Fig.[j]a) , A is equal to the momentum of the undetected 
hadron, p u , while for the radiative process (Fig|lf3,c), it 
is a sum of the undetected hadron and bremsstrahlung 
photon momenta, A = p u + k. 
The phase space is calculated as 



and 



^L^!LS(A- Pu ) = ^dn h = 
2E u 2E h V ' 4W 



8W 2 



(7) 



(8) 



To calculate the matrix element squared, one needs to 
contract the leptonic and hadronic tensors, 



e 4 2e 4 ^ 



(9) 



We consider the longitudinally polarized lepton beam. In 
this case the leptonic tensor reads: 



In the Born approximation, the second term can be 
dropped. However, it gives a non-zero contribution to 
RC. 

For the hadronic tensor, we use a general covariant 
form 

W u.v = -gt*»Hi + P« P U H 2 + pIpIHz (12) 
+ (ra+p£p")W4 + (pKp"-«)W6, 

where the tilde for an arbitrary 4-vector a M denotes the 



substitution a M 



2d nt>- 



g M (to ensure electromagnetic 



gauge invariance). The first four structure functions 
have tensor coefficients symmetric over Lorentz indices, 
but the last one is antisymmetric, it contributes to the 
polarization dependent part of the cross section. The 
quantities 0\ have the following form, 



(13) 



; (SX-M 2 Q 2 ), 

^(SV^ + XV^-S^Q 2 ), 
65 = - 2e ( k i,k 2 ,P,Ph), 
where the operator in the last line is defined as 

In the chosen cm. frame it is equal to 

e(fci, k 2 ,p,p h ) = WP q phPi sin 0\ sin Oh sin <j> h 



(14) 



AW 



(15) 



As a result, we obtain for the Born cross section, (00 = 
da/dW 2 dQ 2 dfl h ) 



32ir 2 S 2 W 2 Q i 1 



(16) 



The expression for the Born cross section in the form 
( ]l6| ) is convenient for our further calculation. It is equiva- 
lent to the well known formula in terms of photoabsorp- 
tion cross sections (or response functions), where each 
term corresponds to certain polarization states of the vir- 
tual photon fl27|, 



L% = ^Tr(fc 2 + m)^(h + m)(l + i^) lv . (10) 

Here the lepton polarization vector is kept in a general 
form. If the lepton is longitudinally polarized and its 
hclicity is positive, then the vector £ can be expressed as 



1 

An 



where 



-co 



or + e &l + e cos 2<fih gtt 



+ y/e(l + e)/2 cos <f> h a LT 
+h ey /e(l-e)/2 eanfa a' LT ,(17) 



1 



5' 



/AT to 



k\ — 2m p) . 



(11) 



W 2 - M 2 
'~2S 2 Q 2 (1- 



(18) 
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h e is the longitudinal polarization degree of the incoming 
leptons, and e |2?J describes the virtual photon polariza- 
tion. 

The following expressions relate the structure functions 
( |l2| ) and Born photoabsorption cross sections ( |l7| ) , 

C(ctt — <jtt), 
2C ( 

2Q 2 ((7 T - (TTT + OX), 



Hi(W 2 ,Q 2 ,t) 

n 2 (w 2 ,Q 2 ,t) 



n 3 {w 2 ,Q 2 ,t) 
n 4 (w 2 ,Q 2 ,t) 
n 5 (w 2 ,Q 2 ,t) 



a 9 . 

—TQa lt + T 2 (Jtt 
2CX q 

-(TTT, 



c_ 
c_ 



{2Tott — Q&lt) 



(19) 



where 



C = 



T 



16ir 2 (W 2 - M 2 )W 2 



A/ = Q 2 S, 



StS x tq 



M 2 t\ 



1 \ 

h A q, 



(20) 



and Q = y^Q 2 . Note that both structure functions Hi, 
as well as cross sections er's, are functions of three in- 
dependent invariant variables, which usually are chosen 
as Q 2 , W 2 and t. Therefore their transformation coeffi- 
cients depend only on these variables but not on <ph ■ For 
the radiative case the variables will be defined in the next 
section (see Eq.(pfl)). For the Born case, they are taken 



as: W 2 , Q 2 and t = t 



v 2 ° 



V° - Q 2 



III. EXACT FORMULAS FOR RADIATIVE 
CORRECTION 



The cross section of radiative process is given by 
M 2 dk 2 dk dp u dph 



da 



2S{2n) s 2E 2 2u 2E U 2E h 
(4ira) 3 dQ 2 dW 2 dn h 



5(A 



2(4ir) 7 S 2 W 2 



dfli-dv- 



PwQ A 



where Qk(h) stands for the solid angle of the 
bremsstrahlung photon (detected hadron). 

Let us first consider the phase space for the radiative 
process. Integrating over the 3-momentum of unobserved 
hadron and using Dirac <5-function to eliminate integra- 
tion over the hadron energy, we have 



dk dp u dp h 
2^2K2El 



5(A-k-p u ) = 



1 



= — / dtt h dQk 




(22) 



where 



fw = W - E h + p h (cos 9 h cos 9, 



+ 9 h sin 9 k cos{4> h - <f> k )). (23) 

Here we introduce a quantity v that describes the missing 
mass (or inelasticity) due to emission of a bremsstrahlung 
photon, v = A 2 — rn^. Note that v = for the Born 
process (Fig|l]a) , as well as in the soft -photon limit (k — 
0) . As can be seen from ( p| ) , both the photon and hadron 
energies are now functions of v. It is related to another 
quantity, fw, as follows: 



v = W 2 + ml 



2WE h = 2ujf v 



(24) 



The largest value of inelasticity allowed by kinematics 
(v m ) corresponds to the threshold of electroproduction. 
It is therefore defined from the relation Eh = mh, yield- 
ing 



(W - m h ) 2 - ml 



(25) 



Note that v m is always smaller for the heavier hardon, 
namely, the nucleon, detected in the final state. It does 
not depend on photon angles, therefore the integration 
region is a rectangle. 

The maximum inelasticity v m is an important quan- 
tity for the RC calculation. All kinematic cuts made by 
experimentalists in data analysis influence RC. It is of- 
ten possible to reduce all these cuts to one effective cut 
on the inelasticity (or missing mass) , v cu t ■ In which case 
Vcut should replace v m as the upper limit of integration 
in Eq. (f22|) , and thus RC can be calculated within the 
cuts using the obtained formulas. If no cuts are applied, 
the maximum value of v equals v m , as given by energy- 
momentum conservation (p5[). 

Now we can fix the kinematics of the radiative process. 
We have to express all scalar products and kinematic 
variables in terms of seven variables: Four variables that 
define the differential cross section and three integration 
variables. First, note that all definitions of leptonic vari- 
ables given in (||) hold in this case. Hadron and real 
photon energies are defined by (24). Hadron momentum 



(21) and scalar products are given as 



Ph = 
V 1>2 = 

St = 



2W ' 



4m? IF 2 , 



\ w = (W -m 2 h - m 2 u - v) - ^,„ h 

= 2(E h2 E h - Pi, 2 p h (cos9 h cos9 1 ^ H 
sin 9 h sin6>i :2 cos 4> h )), 



2pp h = W 2 +m 2 h - ml + V 2 

= S x +t + M 2 - 77 



Vi 



t = V 2 -V 1 -Q 2 



(26) 



As in (||), the last expression is obtained from conser- 
vation laws. All hadron kinematic variables depend on 
inelasticity v but not on photon angles fifc. Scalar prod- 
ucts containing the photon 4-momentum read 



2pk = R w (l — t) = 2lu(E p — p p cos 9k), 



(27) 
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2fci >2 fc = Ki,2 = 2uj(E lj2 ~Pl,2(cOs9k COS 01,2 + 

sin sin 2 cos </>fc), 
2p/jfc = ii = 2uj(E h - ph(cos 9 k cos 6^ + 

sin6> fc sin 6^ cos(</>fc - 4>h))- 

Here we use also invariant variables i?™ = 2k(p + q) and 
t = 2kq/R w ; R w = f w v/W. 

The leptonic tensor of the radiative process has a more 
complicated form. It can be written as 

l>% = ^Tr(/j 2 + m)r M „(fci + m)(l + ijsiftav, 







\l» - 


Inkjet 


la 






2kk x 


2fcfc 2 






)7a< - 








kk\ j 


2kk x 


2kk 2 



(28) 



After contraction of the radiative leptonic tensor, we ob- 
tain 



M 2 R 



2e 6 



(29) 



R w was extracted explicitly in order to cancel v coming 
from Jacobian (p2|). Arguments of Hi can be expressed 
as 



Q 



(q-k) 2 = Q 2 + R w T, 



W 2 = (p + q-kf = W 2 



R-ii 



t = (q - k- p h ) =t - R w (t - fjt). 



(30) 



(31) 



It is well known that the cross section of the radiative 
process is infrared-divergent, which requires careful con- 
sideration in order to cancel in the difference. The pro- 
cedure will be discussed in the next section. Here we can 
extract the infrared convergent terms in 9i in separate 
pieces: 



-5~F IR 9f 



(32) 



The quantities Of are defined by expressions for 0® with 
a reservation that the hadronic quantities Vy , V 2 , St and 
vector ph itself have to be calculated for the radiative 
kinematics (26). This term originates from the first terms 
of r MQ and T av in definition of ( p8[ ) . 

The explicit form of finite parts of these functions Of 
are given in Appendix. 

Finally, the cross section of the radiative process (or = 
d(jR/dQ 2 dW 2 d£l} l ) is given as 



or = 



2 9 7T 4 S 2 W/ 4 



dQhdv 




(33) 



where / = fw/W. 

As a cross-check, we consider the soft-photon limit 
<-"-Vnm < R/2M < Lu max <C all energies and masses. In 
this case only the first term in (32) survives. 

2a f Q 2 \ w max 
or = — log —j - 1 log cr - 

7T V m J Wmin 



(34) 



Integration over photonic angles is performed analyti- 
cally, 

/ dn k F IR = -2(l m - 1), l m = log (35) 



IV. INFRARED DIVERGENCE 

As was already mentioned, this cross section contains 
an infrared divergence. Therefore, in order to compute 
RC, we have to use some regularization method first. We 
use the method of Bardin and Shumeiko for a covari- 
ant treatment of an infrared divergence problem. Basi- 
cally, we follow the original papers devoted to this topic, 
namely, Ref. || for e/i elastic scattering and Ref. U 
where exclusive electroproduction was considered. One 
can find a good and detailed review in Ref. [fL0| . 

Following the rules of dimensional regularization, we 
apply an identity transformation to the radiative cross 
section assuming that we deal with n-dimensional space: 



or = or - oir + oir = a F + oir. 



(36) 



The purpose of the procedure is to separate the cross sec- 
tion into two pieces. The first term doF is complicated 
but infrared free. Infrared divergence is contained in the 
separate term dajR, which has a quite simple structure 
and can be analytically calculated within n-dimensional 
space. In principle, there is some arbitrariness in choos- 
ing the form of the subtracted term. Actually only the 
asymptotical form for R w (or w) — » is fixed. Another 
limitation comes from the theorem about a possibility to 
switch the order of integration and the limit n — > 4. It 
means that we have to provide uniform convergence of 
dap in the limit. Practically, it is ensured if the sub- 
tracted term has a structure F/w, and the difference is 
(F(w) - F(0))/w. 

We define the subtracted part of the radiative cross 
section as follows: Only the first term from Eq. ( |3^ ) 
r.h.s. is kept. It gives the required 1/R W behavior of the 
radiative cross section. Everywhere else, except for the 
5-function, we assume v = 0. This allows to factorize the 
Born cross section when calculating the correction: 



a f f dk 

oir = -o -2 dv —F IR S((A - k) 2 - m 2 J. (37) 

In the cross section of we can now remove the regu- 
larization. The result is 



OF 



(38) 



E 



Q 



For calculation purposes, we split the integration re- 
gion into two parts separated by the infinitesimal value 
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of inelasticity v: 



a rF > a . 
<7ir = -Or co = -{os + o h )vo 

7T 



(39) 



with 

Ss 

5h 



— / dv 

7T 



dv 



F IR S((A-k) 2 
F IR S((A-k) 2 -ml). 



(27T^)"- 4 fc 
d 3 k 

k 



(40) 



The term S S R corresponds to the soft photon contribu- 
tion, while the term Sjj* is caused only by hard photons 
and therefore it does not contain the infrared singularity. 
Again, in the second contribution the regularization can 
be removed. 

We keep integration over the 3-momentum of the radi- 
ated photon in a covariant form. It makes it possible to 
calculate the integral in any frame. We choose the frame 
where A = (so-called R-frame). Integration over in- 
elasticity is external, therefore v is fixed for the integral. 

Calculation of the hard-photon contribution is 
straightforward, 



iff 



(l m - 1) log 



(41) 



For Ss, we follow Ref. || generalizing these calcula- 
tions to the exclusive electroproduction case. 

Using the spherical n-dimensional frame, we have H, 



Sy = - 

IT 



sin 



dv 



2w n/2-l 



(2 7 r At0 )"- 4 r(| - 1) 



dkok^ 



n— 3 



•9d0 



fci k 2 
2k x k ~ 2k 2 k 



S(v - 2k m u )(42) 



It should be noted that the same polar angles are used in 
Eq. (^) for all three terms. That is why we can rotate 
the coordinate system for all the terms and choose the 
OZ axis along k\ t 2,ai respectively. After straightforward 
integration over the polar angle, we obtain Eq. (43). 

The next two steps involve integration over fco using 
the 5-function and over v 



dv I dfcofcg 5 S(v — 2kom u ) — 



1 



n — 4 V 2m u 



n—4 



(46) 

Since the pole (n — 4) 1 is extracted, we can use expan- 
sion in series of n — 4, 



1 



n—4 



-» PlR + l0g ■ 



(tt(1 - xi 2 )) 71 / 2 - 2 



(47) 



with £ = cosQ and standard Pi R defined as 



p, 



i 



i 



IR 



71-4 



:1e 



In 



1 



20?' 



(48) 



"/e being the Euler constant. After then we reduce the 
expression for Ss to the form 



dad(, P IR + log + log Vl-C 2 F(a, £), 

2^i m u 

(49) 

which includes only standard integration that can be 
done using, for instance, tables from Appendix D of 
Ref.jicj]. Integration is straightforward, but it contains 
an integral usually associated with a so-called S$ func- 
tion in the Bardin and Shumeiko approach. Here we skip 
the discussion about properties of the function and give 
the result in the ultra-relativistic approximation which 
we use for the entire calculation: 



F{a,cos{0)) 



k\ k 2 
2k x k ~ 2k 2 k 



4k 2 



rn 



E 2 R (1 - /3 2 cos( 



daQ 2 



S 2 (!-/?« cos ( 



01 COS0) 2 

(43) 



s * = T 



da 



Here Feynman parameters are introduced in order to join 
two denominators. Thus we define the new vector 



k a = aki + (1 — a)k 2 



(44) 



and quantities /3i,2,a &s a ratio of energy to momentum 
in the R-frame. In terms of invariants, we have 




(aS' + (1 - a)X' a y 



(45) 



/3 Q (m 2 + a(l-a)Q 2 ) 
1, 



log 



-Li a 1 - 



Q 2 m 2 u 

S'qXq 



(50) 



Finally, we have for Ss 



5 S = 2 P m + \og 



(J.M 



(lr, 



1) 



log 



S'X' 
m 2 M 2 



+ S^. (51) 



The infrared-divergent terms Pir, as well as the pa- 
rameters /i and v are completely canceled in the sum 
Ss + °~h with Sy which is a contribution of the vertex 
function (Fig.llld): 



-2(P IR +log-)(l m -l)- 



1 2 3 

2 m+ 2 71 



(52) 
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For this sum, we have 
a 

— {os + oh 

where 



where z x = 1 - v/S', z 2 = X'/(X' + v) and 



S v ) = 5, 



a a ( 3 ; 

OVR — — TT'r. 

7T \ Z 



in/ 



(53) 



1 - 



Q 2 M 2 



~6 



Jin/ 



(im-l)ln 



SqXq 



Higher-order corrections can be partially taken into ac- 
count using a special procedure of exponentiation of mul- 
tiple soft photon radiation. There is an uncertainty about 
which part of 5vr has to be exponentiated. Within the 
approach |loj , we have to change 1 + 5j„/ to exp Si n f . 

Collecting all the terms, we obtain for the case of me- 
son electroproduction: 



obs = o- Q e Smf (l + S VR + 5 vac ) + a F . 



(55) 



Here the corrections 5i n f and 5 vac come from radiation of 
soft photons and effects of vacuum polarization, the cor- 
rection Svr is an infrared-free sum of factorized parts of 
real and virtual photon radiation, and op is an infrared- 
free contribution from the bremsstrahlung process. 



W 2 

w 2 



wi 



x o W 2 

A W s,x 



ziS — X — z\i 

S _E_91. 

z-i z 2 



M 



r'2 



M 



(59) 



(54) 



and s x are calculated in accordance with Eq. (Q) 
using W 2 — W 2 X . 

These leading log formulas were extracted from our ex- 
act expressions given in the previous section. The result 
can also be obtained using standard leading log tech- 
niques p2| from our expression (pl|). The leading term 
comes from angular integration of denominators like kki 
and kk 2 . They can be extracted at the level of leptonic 
tensor 



T R T 



1 



1 



1 + Z 



1 



l — 2ii kk\ z 2 {l - z 2 ) kk 2 



(60) 



Only these scalar products in the denominators are sub- 
ject to angular integration. For our exclusive process, 
integration yields: 



V. LEADING LOG APPROXIMATION 

In this section we extract the leading log contribution 
from formulas obtained in the previous section. After 
then we show that the result coincides with what we ob- 
tain from a generalized leading log approximation. 

The method of extraction of leading log (or peaking) 
contributions was first suggested in Ref. 0] (see also 
papers 0). Formally, we have to calculate residues 
of the terms l/z\ and l/z 2 . The corresponding poles 
appear in r for t = t s = —Q 2 /S and r = t x = Q 2 /X. As 
a result, we have (apart from the factorizable correction) 
two contributions of hard radiation, 



cll = — (in 
Ztt 



1) 3 + 21og 



S'X' 



o~a + o~s + o~x 



(56) 



Technically, the contributions can be obtained as 



-8irl r , 
-8ttL 



W 2 



l + z 



S-Q 2 z 1 (l-z 1 
W 2 



1 q b 

i i 



1 + 4 aB 



X + Q 2 I ~ z 2 1 ■ 
For the two peaks the contributions are, respectively: 



dv ( 1 + z\ 2<j Q 
C s a s — 



(57) 



S' V 1 - zx 



z 2 dv ( 1 + z, 



X' 



l-*2 



2 n 



1-ZiJ' 

l-*2 



(58) 



dp u dk 
2e u 2lu 



8(A 



1 



Pu 



fcfel, 



ttI 



(61) 



Now, using (|2^), (|o|) and (|6l|), we obtain the following 
leading-log cross section for the radiative process: 



d 6 cr 



dT b 



al m [ 1 l + z 2 d 6 a (z 1 k 1 ,k 2 ) 
' 2ir [S'l-zt df e 

1 1 + z^d 6 a a (k 1 ,k 2 /z 2 ) 



X'l-Z 2 



dt e 



d 3 k 2 d 3 p h 
2e 2 2e h ' 



(62) 



While deriving the above results, (£7],[58|), we take into 
account relations between z\^ v and S' (X 1 ), which fol- 
low from the constraint (A 2 — (1 — zi)fci) 2 
(A 2 - (1 - l/z 2 )k 2 ) 2 - to 2 = 0. They are 



(l-zi)S', 



1 



Z 2 



X'. 



Z2 



(63) 



VI. 



EXCLUSIVE AND INCLUSIVE RADIATIVE 
CORRECTION 



An important consistency test is to show that inclusive 
RC can be obtained by integration over hadronic angles. 
It should be noted that this is not trivial because the 
hadronic angles of radiative and Born cross sections are 
defined even in different frames. For definition of the 
Born angles, we use the center-of-mass frame, while for 
angles of the radiative process we have to use another 
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frame defined by vectors p and q — k. As a result, it leads 
to quite complicated kinematic relations between these 
angles. 



We start with Eq.(56). After integration over hadronic 
angles, we have to obtain the leading log cross section 
(cr™ = da / dW 2 dQ 2 ) for inclusive case. The inclusive 
formula can be found in fiefs. |3lj ], for example. For 
double differential cross sections in Q 2 and W 2 , it reads 



■mc 



with 



2?r 



{l m - 1) ( (3 + 2 log -He-) < lc + 4" c + 



^1 „mc 

— CTn 



-21 



'OS' 



2a l nc 



(64) 



(65) 



d Zl dQ, s h 



W 2 



S' 



Xw Ziy—D 

w 2 s' v^7V=25; 



The final formulas for the transformation 

v m 1 



(71) 



J dn h J dv = J dzi J dVLl ^j- (72) 



are obtained if we use the equality that can be checked 
directly: 



Ds _ 2 

D 



Ac 



1 



(73) 



2-\ 

2tt 



dz 2 1 + z? 



2 _mc 
On 



'OX 



2a l nc 
I-22' 



(66) 



where z™ = A/u 2 and z™ = U\/S. 

Let us consider integration over fi^ of radiative cross 
sections (i.e., the first terms on the left-hand-sides in 
(|65|)). The simplest way to relate angles in different 
frames is to express them in terms of kinematic invari- 
ants. We use the following relation, 



W 2 



dV?dV 2 ° 



(67) 



where D(Vt\ V%) — [e(fci, fc 2 , p, Ph)] 2 is Gramm determi- 
nant (see (jlj) and Ref . |p3| ) . For the cross sections a % g c 
and u % x c the corresponding expressions are 



dflh = 



W 2 



z 1 dV 1 dV 2 



^X^ s ^-D{z 1 V l ,V 2 ) 



and 



dSl h = 



W 2 



dV 1 dV 2 



(68) 



(69) 



^%^Z 2 y/-D(V 1 ,V 2 /z 2 ) 

The expression for radiative kinematics is 

16D(Vx,V 2 ) = [u 1 V 2 +u 2 V l -{W 2 +ml-ml-v)Q 2 ) 2 
-4(SA-Af 2 Q 2 )(UiU 2 -m£Q 2 ). (70) 

The Born case is reached in the limit v — > 0. Two com- 
ments are in order before explicit integration. The inte- 
gration area in variables V\ and V 2 is defined by equa- 
tion D(Vi,V 2 ) — that produces an ellipse. For radia- 
tive cross sections <t™ c and tr™ c this area is Z\ (or z 2 )- 
dependent. 

A simple integration procedure follows: 



dQh J dv 





dv / dilh 



dV x dV 2 




Using Eq. (|72|), one can see that the exclusive radiative 
cross section transforms to the inclusive one after inte- 
gration. 

Let us show how to obtain a factorized part of the 
inclusive correction from the exclusive result. The sub- 
tracted part of Eqs.(^) can be rewritten as 



dv 



S'(l-zi) 



dz\ 
1-zi 

dz\ 
l-zi 



S'(v=0) 



S'(v=v m ) 

- log 



S' W 



After similar calculation for ax , we can see that inclu- 
sive RC is reproduced exactly. 



VII. NUMERICAL ANALYSIS 

Based on the derived analytical formulas, a 
Fortran code EXCLURAD was developed. This 
code computes RC to the four-fold cross section 
(d A a/dWdQ 2 d cos ddtfih) and to polarization beam 
asymmetry for processes (Q). Both exact and leading log 
formulas obtained in the previous sections are included. 
Any value of the inelasticity cut can be optionally 
chosen. 

Once the computational algorithm for RC is estab- 
lished, there are two possible ways of incorporating RC 
into analysis of experimental data on electroproduction: 

a) Iteration procedure similar to the one implemented 
in the RC code POLRAD for inclusive reactions or 

b) Using realistic models for the structure functions 
dl9| ) of coincidence electroproduction. 

Although the former choice seems attractive due to its 
model independence, it requires, however, full and pre- 
cise experimental mapping of the structure functions ( 1^) 
in the entire range of the kinematic variables needed to 
compute the radiative process (p3|). Such a procedure is 
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0.96 
0.94 
0.92 

0.9 : 
0.88 
0.86 
0.84 ~ 7 
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FIG. 3: VK-dependence of RC to the cross section of neu- 
tral pion production. The models used are MAID2000 (solid 
curve), MAID '98 Q (dashed curve) and AO §j| (dotted 
curve) . 
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1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 

W,GeV 



FIG. 4: VK-dependence of RC to cross section of charged 
pion production. Kinematics and notation are as in Fig.H. 



much more challenging and less efficient than for inclu- 
sive processes. Given that such mapping is not available 
at this time, the only choice left is (b). In this way, 
RC is applied to the model calculations, and then model 
parameters are fixed from available experimental data. 
Thus, RC appears as a necessary intermediate step in 
extracting model parameters from measurements. 

We use the following models for calculation of structure 



FIG. 5: RC to the cross section as a function of cos 9. 



functions ( |l9|) : 

MAID |34| . In this model, baryon resonances are 
described using Breit-Wigner forms, while background 
contributions are described using standard Born terms, 
mixed pseudovector-pseudoscalar 7rNN coupling and t- 
channel vector meson exchange. The final amplitudes are 
constrained by unitarity and gauge invariance. We use 
two versions of MAID: the earlier one, denoted MAID'98 
and the most recent, quoted here as MAID2000. 

AO (Amplitude and Observables) [j35| . The am- 
plitudes are parameterized as follows: S-channel reso- 
nances are parameterized with relativistic Breit-Wigner 
forms with momentum-dependent widths. This part of 
the amplitude is complex. In addition, the s-channcl 
and u-channel pion Born terms are included. These are 
real numbers. Additional real background amplitudes are 
used with an energy dependence that has correct thresh- 
old behavior. 

Another model that can be included into the code EX- 
CLURAD in a simple and straightforward way is the dy- 
namical model of Sato and Lee In this model, the 
off-shell non-resonant contributions to j*p — > A + (1232) 
were calculated directly by applying reaction theory 
within the Hamiltonian formulation underlying 'bare' 
photo-coupling form factors. It should be noted that 
the Sato-Lee model does not include the contributions 
from higher resonances. 

Since the computational algorithm for RC does not 
depend on a particular choice of hadronic structure func- 
tions, addition of any other model does not constitute a 
problem. 

With the models at hand, we can now evaluate nu- 
merically the magnitude of RC as a function of various 
kinematics variables. Let us define the RC factor as fol- 
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FIG. 6: Dependence of RC to the cross section on the az- 
imuthal angle cj>. 



lows, 



5 = 



°~obs 
CTO 



(75) 



For all the following plots, the electron beam energy is 
Ebeam= 1-645 GeV and no cuts on inelasticity were used, 
except for Fig|| We choose representative kinematics for 
current experiments with the CLAS detector at Jefferson 
Lab Hall B. The RC factors to the pion electroproduction 
cross sections are presented in Figs.|| and |asa function 
of W. In this region, the characteristic features of the 
cross section vs. W are strong A(1232) and Sn(1535) 
resonance peaks. One can see that RC is negative at the 
resonance peaks and positive in the dip regions between 
the peaks. Thus it tends to weaken the stength of the 
peaks and fill the dip regions above the given resonances. 
Also shown in Fig. [| is the model dependence effect 
which appears to be noticeable at higher W away from 
the resonance peaks. At higher values of W, the magni- 
tude of RC is larger in the charged pion production case. 
This is mainly due to the wider range of inelasticity v 
( p4| ) associated with detection of a lighter hadron (i.e., 
the pion). 

Before discussing the angular dependence of RC, let 
us first comment on the definitions of hadronic angles 
in order to avoid possible ambiguities. The c.m.s. an- 
gles 9 and <f> are between the direction of momentum lost 
by electrons (<f = k\ — ^2) and momentum of the final 
pion ph, provided that the pion is detected. For the neu- 
tral pion production case, the convention is to also use 
the direction of pion momentum reconstructed from kine- 
matics. In order to follow the convention for the neutral 
pion production, we define 9 and 4> with respect to the 
direction of —pn opposite to the final proton momentum, 



FIG. 7: VK-dependence of the beam polarization asymmetry 
in neutral pion production. The solid (dashed) curve denote 
the asymmetry with (without) RC. MAID2000 was used to 
compute the structure functions. 



while keeping in mind that it is different from the final 
pion momentum for the radiative processes Figjl|b,c. 

The angular dependence of RC is shown in Figs.|| and 
^, where it is plotted as a function of cos# and <fi, respec- 
tively. The kinematics corresponds to the A(1232)-peak, 
where RC leads to suppression of the cross section. It 
can be seen from Fig.|5|that, approaching the forward di- 
rection, at the given value of <j), the magnitude of RC as a 
function of cos 9 for the charged-pion channel smoothly 
increases from about 3% to 12%, while for the neutral 
pions is decreases from 10% to 5%. A common feature 
for both the channels is that RC is larger in magnitude 
for the parallel kinematics, when the detected hadron 
moves along the transferred momentum. The RC factor 
varies as a function of azimuthal angle 0, as well (FigJ^). 
Dependence of RC on the angle cf> has important implica- 
tions for the super-Rosenbluth separation of electropro- 
duction structure functions, and the 9 dependence would 
affect the partial-wave analysis, resulting in corrections 
to electroexcitation parameters of baryon resonances. 

The beam polarization asymmetry is plotted in Figs. 
and |^ for the neutral and charged channels, respectively. 
One can see that for the asymmetry, RC changes from 
enhancement to suppression when passing across the A 
and Sn resonance regions. The RC factor is most sub- 
stantial in the dip regions between resonances. 

Figure |p demonstrates the dependence on inelas- 
ticity cut v cu t (24). It can be seen that for the 
smaller values of the cut, resulting in selection of 



softer bremsstrahlung photons, RC to both polarization- 
dependent and polarization-independent parts of the 
cross section is almost the same, resulting in a small cor- 
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FIG. 8: The beam polarization asymmetry in charged pion 
production as a function of W with (solid curve) and without 
(dashed curve) RC. Notation is as in Fig.[|. 



rection to polarization asymmetry. As v cu t increases, the 
correction to the asymmetry also increases due to the 
hard-photon emission coming into play. 

Any experimental spectrum for electroproduction con- 
tains a radiative tail that, due to the finite energy resolu- 
tion, cannot be experimentally separated from Born con- 
tribution. In practice, detected events are included into 
data analysis up to a certain (cut) limit and are inter- 
preted as radiative events. The contribution of radiative 
events should be calculated theoretically using the same 
inelasticity cut (v cu t). The resulting contribution of the 
radiative events is subtracted from the integrated exper- 
imental spectrum in order to obtain the value of Born 
cross section that, naturally, should be cut-independent. 

The Fortran code EXCLURAD can be downloaded 



from http: / /www .jlab.org/RC or obtained directly from 
the authors. 



VIII. DISCUSSION AND CONCLUSION 

In this paper, we obtain explicit formulas for the 
lowest-order QED radiative correction to cross section 
and polarization beam asymmetry in the exclusive pion 
electroproduction. Analytic formulas are tested in sev- 
eral ways. Apart from traditional cross-checks like soft- 
photon and leading-log limits, it is found that integra- 
tion with respect to the hadronic angles reproduces the 
inclusive radiative correction. 

A Fortran code EXCLURAD is developed on the basis 
of the analytic formulas. Numerical analysis carried out 
for JLab kinematic conditions shows that: 



FIG. 9: Dependence of RC to cross section and polarization 
beam asymmetry on the inelasticity cut v for neutral pion 
production. The quantity 8 unp is RC to the unpolarized part 
of the cross section, S po i is RC to the polarized part of the 
cross section, and 5a = S unp /5 po i, i.e., it is RC to the beam 
polarization asymmetry. MAID2000 was used for structure 
functions. 



• Radiative correction to the cross section of electro- 
production is very sensitive to the cut on inelastic- 
ity. The harder cut eliminates contributions from 
the higher-energy part of the bremsstrahlung spec- 
trum, thus leading to the smaller magnitude of RC 
to the polarization asymmetry. 

• RC to cross sections can be as high as several tens 
of per cent. 

• RC may vary depending on the chosen model for 
electroproduction structure functions. The pro- 
posed RC procedure may be viewed as a necessary 
intermediate step in interpretation of experimen- 
tal results in terms of model parameters. Iteration 
procedure may be required for the regions where 
model dependence is substantial. 

• RC have a non-trivial angular dependence in cos 9 
and 4>. This is particularly significant, as these are 
often used as input to extract structure functions 
and partial wave amplitudes. 
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APPENDIX A 

In this Appendix, we give formulas for Of ( Eq. (|32| 
In the unpolarized case, 



Q[ = ®i1 + Rw@i3i (Al) 

-AF-2F d ,T 2 , 

-AM 2 F ir t - F d S 2 p T + F 1+ S P S X 
+2F 2 _S P + 2Fj R S x , 

AM 2 F + 2M 2 F d r 2 - F d S x r - F 1+ S p , 

-Am 2 F IR T - F d rV 2 + F 1+ V-V+ 
+2F 2 _ t iV + + 2 f iF IR V_, 

4m 2 F + 2m 2 F d T 2 - F df irV- - F 1+f iV+, 

-2F d S p rV + + F 1+ S P V- + F l+ S x V+ 
+2F 2 _ f iS p + 2F 2 _V+ + 2[mF ir S x 
-AF IR S m T + 2F IR V_, 

AFS m - F d nS x r + 2F d S m r 2 - F d rV- 
-F 1+ ^S P - F 1+ V+. 
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26*22 

2^23 
2^32 

26*33 
26*42 

26*43 



Here V± = V l ±V 2 , 
F d 
Fi+ 



K\K 2 
1 1 

«1 K 2 



(A2) 
(A3) 



F 2 ± 
Fir 



m " 



F, 



2+ 



± 



Q 2 F d 



(A4) 
(A5) 



where K\. 2 are defined in Eq. (Efl). 



7.V2 



7.13 



Al^F^ + 2F 22 E l + F d (E l2 r - Q 2 (E 1 + E 2 )) 

+ ^(E 1 S - E 12 t + E 12 - E 2 S)j, 

2{f 1+ {E 2 - Ei) - F d {E t + E 2 )t 
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F tl E 2 {r-l) ), 



(A6) 



E12 



E, = 



e(h,k 2 ,pi,Ph) 



1 



■(V 1 X-Q*S m + SV 2 ) 
-i 1/2 

+A{SX-M 2 Q 2 ){Q 2 m 2 h -ViV 2 ) 



e(k,ki,pi,ph) = j 



-{RMt - 1) - i^/zS + S mKl ) 2 

-,1/2 



+4(kiAP + R w (t - l)S)(-Kim% + ViRun) 



E 2 = 



e(k,k 2 ,pi,Ph) = j 



-CR»Va(T - 1) - Rw^X + S m n 2 ) 2 
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+4(/« 2 M 2 +i? w (T - l)X)(-K 2 m£ + V 2 RwfJ,) 
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